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ABSTRACT 

We continue to explore the validity of the reflected shock structure (RSS) picture in 
SNR f 987A that was proposed in our previous analyses of the X-ray emission from 
this object. We used an improved version of our RSS model in a global analysis of 14 
CCD spectra from the monitoring program with Chandra. In the framework of the 
RSS picture, we are able to match both the expansion velocity curve deduced from 
the analysis of the X-ray images and light curve. Using a simplified analysis, we also 
show that the X-rays and the non-thermal radio emission may originate from the same 
shock structure (the blast wave). We believe that using the RSS model in the analysis 
of grating data from the Chandra monitoring program of SNR 1987A that cover a 
long enough time interval, will allow us to build a more realistic physical picture and 
model of SNR 1987A. 
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1 INTRODUCTION 



> ■ After 



its reappeara nce in 



' Istavelev-Smith et al.l 
rBeucrmann et al.l 



1992I . 



IHii; 



radio (JTurtle et all Il990l : 
1993) and in X-rays 
iGorenstein et al.l 1 19941 : 



(Istavelev-Smith et all Il992l . 1 993, 2007; iBall et all 1200 ll : 
Manchester et all |2002|. 2005: [Caensler et all Il997l. 2 007: 



JHasingcr ct all Il996l ') at about 1000-1200 days after 
the explosion (DAE), SN 1987A entered into the phas e 
of supernova remnant (SNR) formation IjMcCravl 120071 ). 
Observations of this phenomenon from its onset give us a 
chance to study in detail the physics of strong shocks as well 
as to investigate the distribution of the circumstellar matter 
(CSM) around the exploded star. The latter may help us 
solve the mystery of t he origin of the triple - ring system ob- 
served in the optical ([Wampler et al. 1990l: Jakobsen et al.l 
ll99ll : ICrotts fc Heathcotelll99ll : lBurrows et allll995h . With 
this respect. X-ray observations are very important since 
they manifest most of the energetics at present and they 
provide us with direct information about the underlying 
physics of this evolving object. 

The richness and variety of processes involved in the 
birth of SNR 1987A are manifested by observations across 
the entire electromagnetic spectrum. Radio observations 
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Ng et al.ll2008l : JPotter et alll2009l : IZanardo et alll2010l ) tell 
us about the non-thermal (synchrotron) emission mecha- 
nisms and relativistic particle acceleration in strong shocks. 
Near-infrared observations tell us about the distribution 
of dust in the CSM and the physics of its interaction 



with the hot plasma (|Bouchet et al.ll2003 . 2006: 'Pwe k et al 
12008^. O ptical and UV observations (Sonncrborn et al 



^1998; La wrence et al.ll2000l:lPun et alll20o"2l : ISugerman et al 
12002, : iGroningsson et all |2006| . 2008) tell us about the 
physics of slow shocks and the distribution of dense clumps 
of gas. 

The early evolution of SNR 1987A is strongly influ- 
enced by the distribution of circumstellar gas, which is con- 
centrated in its triple-ring system. Although the origin of the 
ring system is not well understood, the colliding wind sce- 
nario that assumes a spherically-symmetric blue supergiant 
(BSG) wind interacting with an earlier emitted asymmetric 
red supergiant (RS G) wind seems the physically most rea- 
sonable explanation (Luo & McCray 1 99ll : IWang fc M azzali 
1992; Morr is fc Podsiadlowski 20 7). It has been explored 
in detail bv JBlondin fc Lundgvistj (|l993l ) who demonstrated 
that such a model is capable of accounting for the observed 
CSM characteristics, provided the RSG wind was highly 
asymmetric with most of its mass confined near the equa- 
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tor. This picture needed an additional ingredient in order 
to ex plain the reappearance o f SNR 1987A in radio and X- 
rays. IChevalier fc Dwarkadaa ([1995!) proposed that an ex- 
tended HII region was formed inside the inner equatorial 
ring that was naturally produced by the photoionization of 
the shocked RSG wind (the inner ring) by the central pro- 
genitor BSG star. Based on this, Borkowski et al. (1997a, 
b) studied the complex hydrodynamics of the interaction of 
the SN ejecta with the GSM in vicinity of the inner ring that 
also predicts a continuous brightening of SNR 1987A in UV 
and X-rays. 

Since its reappearance, SNR 1987A has been brighten- 
ing continuousl y in X-rays with an accelerat ing rate in the 
past few years (|Park et al.ll2005l . 2006, 2007; iRacusin et al] 
|2009| ). The first spectrum of SNR 1987A with high spec- 
tral resolution was obtained with the high energy trans- 
mission grating on Chandra by Mi chael et al.l (|2002f ). The 
subsequent brightening has enabled dispersed spectra with 
good photon statist ics to be obtaine d with both Chandra 
fehekovet al]|2005l . 2006: iDewev et al. 2008; Zhckov et al 



how they might be related to the non-thermal radio emission 
from SNR 1987A and the consistency of the global RSS 
picture. We present our conclusions in section § [T] 



2009) and XMM-Newton l|Haberl et al.l 120061 : iHeng et al 
20081 ). With the superb spatial and spectral resolution of 



the Chandra observatory, the first direct evidence was ob- 
tained that the distribution of the X-ray emitting plasma 
in SNR 1987A is not spherically-symmetric but it is 'flat', 
i.e., it is concentrated in th e vicinity of the e quatorial r ing 
JZhekov et al. 2005 . 2006; Dewey etal .1120081 : [Zhckov ctlo] 
I2OO9I ). These authors also showed that the bulk gas velocities 
deduced from the X-ray spectral lines are too low to account 
for the plasma temperatures inferred from the spectral fit, 
and showed that a simplified reflected-shock structure (RSS) 
model could match the observed spectra successfully. 

Superb sp atial resolution of Chandra allowed 
iBurrows et al.l i|2000l ) to obtain the first image of 
SNR 1987A that showed a ring-like structure in X- 
rays. Subsequently, a monitoring program was established 
with Chandra to track the evolution of the X-ray images 
and GCD spectra of SNR 1987A. In addition to the 
continuing brightening in X-rays, the spectral data showed 
a gradual decrease of the plasma temperature behind the 
fast shocks. From the imaging d ata, we measure d the 
expansion velocity of the remnant (|Park et al.ll2002l . 2004, 
2007). After applying a ph ysically more realisti c method to 
analyze the X-ray images, 'Ra cusin et al.l (!_2009) obtained a 
more accurate expansion velocity curve. It showed clearly 
that the velocity started to decelerate after ~ 6000 DAE, 
about the same t ime that the X-ray light curve turned up 
iPark et al.ll2005l ). 

Here we analyze the evolution of the X-ray spectrum 
of SNR 1987A as observed with the pulse-height spectra 
obtained with the AGIS GCD camera during the monitoring 
program. One of the main goals of our study is to test the 
refiected-shock structure picture that emerged from analysis 
of the Chandra grating data. Since the pulse-height spectra 
have much better time coverage than the grating data, they 
can test whether is possible to match the CGD spectra and 
the expansion velocity curve in the framework of the RSS 
picture. In section §[2] we give a qualitative description of the 
RSS picture for SNR 1987A. In section §[3]we present details 
about the global RSS model. We describe the evolution of 
the GCD spectra and results from the spectral fits in sections 
§[4]and §[5l respectively. In section §[6]we discuss the results. 



2 PHYSICAL PICTURE AND MOTIVATION 

The physical picture that emerges from the analysis of 
SNR 1987A so far is the following. The circumstellar matter 
around the inner equatorial ring likely consists of at least 
two components: (1) a smooth gas component which repre- 
sents the so called HII region (resulting from evaporation of 
the inner ring or being the hot bubble in the interacting stel- 
lar winds scenario); and (2) much denser clumps distributed 
within the smooth component. Although the origin of the 
dense clumps is unclear, it seems realistic to expect that they 
wil l form as a result from vari ous dynamical instabilities. 
As iBlondin fc Lundgvistl l| 19931 ) point out, in the interact- 
ing stellar winds scenario these instabilities are expected to 
produce a very clumpy and nonuniform shell. Thus, during 
the shell evolution dense clumps penetrate into the smooth 
component and form protrusions of various shapes. 

Interaction of the SN ejecta with the smooth component 
gives rise to the standard two-slab structure bounded by 
a forward shock (blast wave compressing the smooth GSM 
component) and a reverse shock (compressing the ejecta ma- 
terial). The slabs are separated by a contact discontinuity. 

When the blast wave overtakes a clump, it is 'lit up', 
namely, a transmitted shock starts to heat up the clump 
gas. At the same time, a reflected shock is sent back that 
shocks the gas behind the blast wave again, further increas- 
ing its temperature and decreasing its bulk velocity. Such 
a reflected-shock structure was invoked to explain the rela- 
tively narrow spectral lines detected by the Chandra grat- 
ing observations of SNR 1987A and the corresponding sim- 
plifled model was successful in describing these results. It 
should be kept in mind that the actual hydrodynamics of 
such an interaction is very complex: the clumps may have 
various shapes; the blast wave may wrap up around a clump 
and result in a more complicated shock structure; the re- 
flected shocks from adjacent clumps may additionally inter- 
act to create a 'turbulent' system of gas shocked multiple 
times. A proper description of such interaction would re- 
quire 3D hydrodynamic calculations. Thus, the RSS model 
might be considered as a very simple representation of the 
actual physical picture. 

On the other hand, it does not seem feasible to attempt 
a complete match of the SNR 1987A observables (X-rays 
spectra, light curve etc.) by a detailed hydrodynamic model 
mostly because there are too many initial (and boundary) 
conditions to explore. Instead, we think it would be instruc- 
tive to use a simplifled model that bears the basic physics of 
this complex interaction and by using such a tool we hope 
to deduce valuable information about the GSM in the vicin- 
ity of the inner ring of SNR 1987A. Results from such an 
analysis of the available (and future) observations can guide 
us towards building a more realistic picture (and model) of 
this exciting phenomenon: the formation and development 
of SNR 1987A. 
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HII region; pRi„g and po are the densities respectively at 
these radii; D and s control the shape of the density profile. 
We assume that the outer supernova envelope has a den- 
sity profile, p oc t~^(r/t)~" , where w = 9 in the case of 
SNR 1987A (e.R., lEastman fc Kirshnerl [l989|: ISuzuki et all 



1 19931 : iBorkowski et all Il997al .b). IChevaUerl (l982l ) derirod 
the solution for the structure that results from interaction 
of such a density profile with the circumstellar medium. 
Guided by his approach (i.e., using dimensional analysis), 
we derive an analytical solution for the radius, r, of the 
structure (the contact discontinuity). In the case of the CSM 
density profile given by eq. |T| we find: 



Figure 1. The analytical solu tion fits to the radial-expansion 
curve from lRacusin et al.l 1 120091 ). The solid line is an overlay of 
a few solutions having tmax from 10,000 to 15,000. The corre- 
sponding density profiles of the HII region are shown in the inset 
figure. 



3 GLOBAL RSS MODEL 

IZhekov et al.l (|2009l ) proposed a reflected-shock struc- 
ture model to interpret the Chandra grating data for 
SNR 1987A taken in 2004-2007. It bears the basic physics 
of interaction of a shock with a dense clump which is a cor- 
nerstone in the physical picture described above. Since the 
CCD monitoring observations provide more extensive time 
coverage of the SNR 1987A evolution, they can be used to 
test details of this physical picture. To do so, we adopt the 
following improvements of the RSS model: 

• At every moment (observation) the entire (integrated) 
X-ray characteristics of SNR 1987A can be represented 
with a RSS structure (blast wave, transmitted and reflected 
shocks), i.e., the blast wave is impinging on some average 
clump that is representative of the clumps at that radius 
from the center of the ring (the mean of the clumps distri- 
bution). 

• The emission measure (EM) of the blast wave is not a 
free parameter from one observation to another but it follows 
the density profile of the smooth component (EM oc n^r^ , 
where n is the density of the smooth component and r is the 
radius of the SNR 1987A at that moment). 

• The SNR 1987A radius is calculated from ID analytical 
solution that takes into account the density profile of the 
smooth component. 

• The electron and ion temperatures are not equal behind 
the shocks (but see below) . 



3.1 Analytical solution 

We anticipate that the density profile of the smooth compo- 
nent will be increasing steeply towards the inner ring radius 
as measured from the Hubble Space Telescope (HST) images. 
Many choices of such a function are possible and the model 
chosen here is: 



T — Tr] 



p = po + [pRtng - po) tank D 



Ro 



Rmng — f J 

where Raing = 0.83 arcsec IjSugerman et al.l [2003) is the 
radius of the equatorial ring; Rq is the inner radius of the 
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where po = po/pm-ng and tmax is the time when the radius 
of the blast wave equals that of the ring. 

We found that this expression for r is in satisfactory 
agreement with the result of a ID hydrodynamic simulations 
for the case of stellar ejecta with n — 9 interacting with CSM 
that has density profile as given in eq. ([T|. This makes us 
confident in using the analytical solution in our model. 

Figure [l] shows that equation ^ gives a good fit t o 
the observed radial-expansion curve ( Racusin et al.ll2009l V 
Note that there is not a unique set of fit parameters {tmax, 
po,D and s) that represents the 'best' fit. Namely, parame- 
ters tmax and Po cross talk and this is illustrated in Fig. [T] 
where the solid curve is in fact an overlay of a few solutions 
having tmax in the range from 10,000 to 15,000 DAE (the 
corresponding density profiles are shown in the inset figure). 
This behavior is well understood since the fit to the radial- 
expansion curve is determined by the relative change of the 
density of the smooth component between the observations. 
It will become possible to discriminate among different cases 
as the blast wave progresses further into the equatorial ring. 



3.2 XSPEC model 

Based on the analyti cal solution, we ext ended the RSS 
model for XSPEC Fl bv IZhekov etakl (|2009l '). The new RSS 
model tracks the radiu s of the blast wave that is inferred 
from the observations (jRacusin et al.ll2009l ). The emission 
measure behind this shock is scaled correspondingly (EM oc 
n^r^) and its ionization time (uet) follows the density pro- 
file of the preshock gas (in fact, its integral over time for a 
particular observation). Also, the blast wave velocity is esti- 
mated for each observation. Bearing in mind that the blast 
wave vel ocity has decreased recently from a few thousands 
km s~^ iRacusin et al.ll2009l ) and the fact that the electron 
temperature of the X-ray emi tting plasma was not greater 
than 3 keV IjPark et al.ll200q ). we conclude that the elec- 
tron and ion temperatures in the plasma behind the blast 
wave have not completely equilibrated (the mean plasma 
temperature kT= 1.4 x [Vshock/WOO km s"^]^ keV for the 



1 Version 11.3.2 of the XSPEC analysis package dArnaudlllOOl) 
is used in this study. 
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SNR 1 987A abundan c es). W e recall that in the original RSS 
model, IZhekov et alj (|2003 ) assumed a complete tempera- 
ture equilibration (IT approximation) and took into account 
the non-equilibrium ionization (NEI) effects explicitly, ffere 
we adopted the same approximation which can be justified 
by the following considerations. 

First, a direct comparison between the X-ray spectra 
from a plane-parallel shock (PPS) model with different elec- 
tron and ion temperatures (2T approximation; vnpshock in 
XSPEC ) and from a IT PPS model {vpshock in XSPEC ) 
shows that the former model can be matched well by the 
latter having a postshock te mperature averaged ov er the 
shock ionization age (see also iBorkowski et alJl200ll ). This 
simply indicates that the NEI effects are more pronounced 
than the effects due to the electron temperature evolution 
downstream from the shock front. Second, we constructed a 
RSS model that correctly treats the temperature equilibra- 
tion behind the shock (similarly as in vnpshock in XSPEC ). 
But we found that its application to the entire set of mon- 
itoring observations is very unstable because the parameter 
space is too complicated. 

Thus, the new RSS model correctly estimates the mean 
plasma temperature behind the blast wave. That result is 
then used to calculate the mean plasma temperature be- 
hind the transmitted sho ck (as well as behin d the reflected 
shock: see Appendix B in lZhekov et al.ll2009l ). The spectral 
fits assume an average electron temperature (Te) typical for 
either of the shocks. From these spectral fits we derive a 
parameter which is the average /3 — T^/T {T is the mean 
plasma temperature) for each shock and its value sets an up- 
per limit to the actual /3 in the shock front. We note that the 
refiected shock is assumed t o have the same ave rage /3 as the 
blast wave. As discussed bv lZhekov et al.l (|2009l ). the plasma 
behind the refiected shock has quite complicated heating 
history. It is so since the refiected shock interacts with the 
plasma behind the blast wave that has already evolved for 
some time. As a result, in the region behind the reflected 
shock we have a mixture of hot plasma with different heat- 
ing histories, resulting in a distribution of temperatures and 
ionization time scales. To simplify the model calculations 
and also due to the technical problems mentioned above, we 
adopted an average jS parameter for the reflected shock that 
is equal to that in the blast wave. 



4 CCD SPECTRA 

The CCD spectra from the monitoring program of 
SNR 1987A are exactly the same (source and background 
sp ectra, redis t ributi on a nd ancillary respo nse functions) as 
in lParket all l|2006l ) and iPark et"al] (|2007l ). Details of data 
reduction and spectral extraction are found therein. In this 
work we analyze 14 data sets, with Chandra Obslds as fol- 
lows: 122, 1967, 1044, 2831, 2832, 3829, 3830, 4614, 4615, 
5579 and 6178, 5580 and 6345, 6668, 6669, 7636 (two of the 
observations were split into two subsets). 



5 SPECTRAL FITS 

As in previous studies of the Chandra spectra of SNR 1987A, 
we assumed that the X-ray absorption and abundances of 



Table 1. The RSS Model Results (CCD Spectra) 



xVdof 


1706/1473 






Nh(1021 cm-2) 


2.23 ±0.05 






H 


1 


1 




He 


2.57 


2.57 




C 


0.09 


0.09 




N (1.63) 


0.32 ±0.04 


0.56 [0.50 - 


■ 0.65] 


O (0.18) 


0.07 ±0.01 


0.081 [0.074 - 0.092] 


Ne(0.29) 


0.22 ±0.01 


0.29 [0.27 - 


■ 0.31] 


Mg(0.32) 


0.17 ±0.01 


0.28 [0.26 - 


■ 0.29] 


Si(0.31) 


0.26 ±0.01 


0.33 [0.32 - 


■ 0.35] 


S (0.36) 


0.36 ±0.02 


0.30 [0.24 - 


■ 0.36] 


Ar 


0.537 


0.537 




Ca 


0.339 


0.339 




Fe(0.22) 


0.10 ±0.01 


0.19 [0.19 - 


■ 0.21] 


Ni 


0.618 


0.618 





Note - Results from the simultaneous fit to 14 CCD spectra 
of SNR 1987A. The uncertainties are la errors from the fit. 
All abundances are expr essed as ratios to their solar values 
IjAnders fc Grevesselll989l ) . For comparison, the inner- ring abun- 
dances of He, C, N, and O JLundavist &: Franssonll996h and those 
of Ne , Mg, Si, S and Fe typical for LMC SNRs JHughes et all 
119981) are given in the first column in parentheses . The Ar, Ca 
and N i abundances are representative for LMC (Ru ssel fc Dopital 
Il992n . The derived abundances (with 90% confidence interval 
in brackets) from t he analysis of the Chandra grating spectra 
llZhekov et al.ll2009l ) are given in the third column. 



the hot plasma do not change with time. Thus, all spectra 
shared the value for Nh and had the same abundances. The 
procedure s for fltting abundan ces are the same as those de- 
scribed in I Michael et al.l ()2002l ). Namely, we only varied the 
abundances of elements having strong emission lines in the 
observed (0.5 - 4 keV) energy range: N, O, Ne, Mg, Si, S and 
Fe. 

To fit the global RSS model to the data we adopt the 
following procedure. We choose one of the density profiles 
in the HII region that provides a good fit to the expansion 
velocity curve (Fig. [T]). For the results presented here we 
used the case with t^ax = 13, 000 DAE and po = 0.0289 
(po = Po/pRing)- Wc held thcsc parameters fixed in the fits. 
Consequently, we know the evolution of the emission mea- 
sure (EM ex v?r^) and ionization time (uet) of the blast 
wave, but in relative units. We derive and their absolute cal- 
ibration from the fits. From the analytical solution we also 
know the velocity of the blast wave which in turn deter- 
mines the mean plasma temperature behind this shock. For 
each observation, we fitted for the average /3 which is a free 
parameter. We derive the mean plasma temperature in the 
tran smitted shock s elf-co nsistently as described in Appendix 
B in lZhekov et al.l (|2009l ) and we only fitted for the average 
electron temperature in this shock. For each observed spec- 
trum we fitted both the emission measure and ionization age 
of the transmitted shock. The RSS model also estimates all 
the parameters (and the spectrum) of the reflected shock in 
a self-consistent way. 

The global RSS model was used to flt 14 CCD spec- 
tra simultaneously. Figure [2] shows that the RSS model 
is very successful in representing the X-ray emission from 
SNR 1987A. Other fit resuhs are given in Table [H Figs. [3] 
and[4l One can see that some shock parameters follow the 
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Figure 2. The background-subtracted spectra of SNR 1987A overlaid with the model fits for the case with tmax = 13, 000 DAE (see 
§|5ll. Each panel is marked by the days after explosion of the corresponding observation. All the spectra are re- binned to have a minimum 
of 50 counts per bin but the observations at 4711 DAE (20 cts/bin) and at 5178 DAE (30 cts/bin). 



same pattern as derived in previous analyses (e.g. lPark et al.l 
120061 ) but their values are derived now in the framework 
of a (highly simplified) physical model. For example, the 
average electron temperature in the blast wave decreases 
gradually, while the mean electron temperature behind the 
transmitted shock increases slightly. These results are con- 
sequences of the increasing preshock density and the corre- 



sponding decreasing shock velocity, which helps the electron- 
ion temperature equilibration. As Figure |4] shows, the blast- 
wave emission measure increases smoothly at an accelerat- 
ing rate. In contrast, the emission measure of the transmit- 
ted and reflected shocks increases suddenly for epochs after 
~ 6200 (6000-6400) DAE in coinciden ce with the steep up- 
turn in the X-ray LC (|Park et al.ll2005h and the decrease in 
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Figure 3. Evolution of various shock parameters as derived from the global RSS model fit. Shown are: the average electron temperatures 
for the blast wave, transmitted and reflected shocks (see text for details). The lower right panel represents the ratio of the average electron 
to the mean postshock temperature for the blast wave (squares) and the transmitted shock (asterisks). For comparison the blast wave 
velocity is denoted with dots and a dotted line, while the velocity of the transmitted shock is depicted with dots and a dashed line. The 
error bars in the temperature plots are la errors from the fit. BW, TS and RS stand respectively for the blast wave, transmitted and 
reflected shock. 



the e xpansion velocity curve of SNR 1987A (JRacusin et al.l 
|2009| ). This behavior is probably a consequence of the fact 
that the number of dense clumps overtaken by the blast wave 
started to increase at those times. Future observations will 
tell us about their distribution deeper in the equatorial ring. 
Also, Figure |4] shows the ionization ages (uet) for each RSS 
component. We note that the high values for the transmitted 
shock in the early observations likely result from numerical 
uncertainties due to the small contribution of this compo- 
nent to the total X-ray emission. Thus, we consider these 
(uet) values less reliable. 

Note that on average the derived abundances have 
smaller value s than the ones de r ived from analysis of the 
grating data l|Dewev et al.ll2008l : IZhekov et aLll2009t ). This 
result is likely an artifact of the diminished resolution in 
the CCD spectra. Figures [5] and [6] show how the RSS 
model matches the X-ray light curve in the 'soft' (0.5-2 
keV) an d 'hard' (3- 1 keV ) en ergy bands . All t he fluxes 
are from iPark et al.l l|2006l ) and iPark et~aLl (|2007l ). We see 
that the emission from the blast wave dominated the soft 
X-ray hght curve until the turning point mentioned above, 
~ 6200 (6000 - 6400) DAE, after which time the emission 
from the transmitted shock became dominated. On the other 
hand, the blast-wave emission continued to dominate the 
hard X-ray hght curve after the turning point, augmented 
by an increasing contribution from reflected shocks. 

Interestingly, the (0.5-2 keV) model flux, extrapolated 



back to ROSAT times, matches those SNR 1987A fluxes 
as well. This extrapolation assumes that at that time the 
entire X-ray emission was formed in the blast wave. There- 
fore, the blast wave flux from the first Chandra observa- 
tion was scaled back with the emission measure of the blast 
wave (EM oc n^r^) which is known from the analytical solu- 
tion. This result gives us more confidence in the global RSS 
model. 

Figure |3] shows that the average /?(/? = Te/T) for 
the blast wave gradually increases as the shock velocity de- 
creases. This behavior is similar to the one for the ratio of 
the electron and ion temperature set at the shock front be- 
ing reversely proportional to the Alfv en Mach number of 
the shock (e.g. iGhavamian et al.ll200ll ). The velocity evolu- 
tion of the average /3 can likely be attributed to this fact as 
well as to the increasing ionization age of the blast wave. 
However, one should keep in mind that the actual physical 
picture is not that simple since any new gas parcel that is 
overtaken by the blast wave does not follow exactly the evo- 
lution of the previous gas parcels and has its own evolution 
instead. 

Note also that the average /3 of the transmitted shock 
follows exactly the same pattern as that for the blast wave. 
We regard this result with caution: it might just be a coin- 
cidence that the values of the two /?'s are about the same 
since although relatively high the velocity of the transmitted 
shock is considerably lower than that of the blast wave. 
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Figure 5. The soft X-ray LC of SNR 1987A from ROSAT and 
Chandra. The data with error bars are from Park et al. (2007). 
Components of the RSS model are marked with; asterisks (the 
blast wave); diamonds (the transmitted shock); squares (the re- 
flected shock). The solid line denotes the extrapolation from the 
first Chandra observation back to the ROSAT times under the 
assumption that all the X-ray emission is due to the blast wave 
alone. 
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Figure 4. Emission measure of the shock components as derived 
from the global RSS model fit. BW, TS and RS stand respectively 
for the blast wave, transmitted and reflected shock. The corre- 
sponding ionization ages are shown in the inset figure of each 
panel. 



Summarizing the basic results, we note that the global 
RSS model provides a simple unified physical picture that 
can match the observed expansion of the X-ray image as 
well as the X-ray light curve and spectral evolution of 
SNR 1987A. We are encouraged to believe that by using 
the RSS model to analyze the data from the ongoing mon- 
itoring of this object with Chandra we may deduce further 
details, not only about CSM, but also about the kinemat- 
ics of the X-ray emitting plasma. Moreover, we are hopeful 
that the information derived from the RSS model analysis 
of the X-rays from SNR 1987A may be used in conjunction 
with observations at other wavelength bands (e.g., radio, 
infrared) to further constrain the model. 



Figure 6. The hard X-ray LC of SNR 1987A from Chandra. 
The data with error bars are from Park et al. (2007). Components 
of the RSS model are marked with: asterisks (the blast wave) and 
squares (the reflected shock). The contribution of the transmitted 
shock is very low and thus it is out of the shown range on the 
y-axis. 



6 DISCUSSION 

6.1 X-rays vs. Radio Emission 

X-rays from SNR 1987A are of thermal origin and they 
provide us with direct information about the strong shocks 
in this object. On the other hand, the radio emission 
from SNR 1987A is non-th ermal synchrotron radiation (e.g. 
IStavelev-Smith et al.ll 19921 ) . which requires two ingredients: 
relativistic electrons and a magnetic field. The relativistic 
electrons are assumed to be born in strong shocks through 
the mechanism of difi^usive shock acceleration (DSA). Thus, 
the non-thermal (NT) radio emission can provide us with 
additional details about the physical picture in SNR 1987A. 
An interesting question is whether the X-rays and NT 
radio emission originate from the same shocks. In what fol- 
lows, we assume that this is indeed the case, specifically, that 
the radio emission comes from the blast wave in SNR 1987A. 
To explore this possibility, we adopt an approach similar to 
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the one we took in the analysis of the CCD X-ray spectra 
above. Namely, we adopt a simple model for the NT radio 
emission and derive some parameters of the relativistic elec- 
trons involved in this process. Because it is not possible to 
measure the magnetic field directly, we assume that it is 
proportional to the density of the shocked plasma (frozen-in 
conditions) . 

The light curves (LC) of SNR 1987A at fo ur radio fre- 
quencies were taken from I Zanardo et al.l (|201Cl ). Under the 
assumption that this is synchrotron emission from relativis- 
tic electrons having a power-law distribution in energy, the 
LCs were fitted simultaneously with a power-law function 
whose normalization and power-law index are function of 
time: 



F.{t) ^ Fo{t) u^^n 



GHz 



(3) 



where i^Qjj^ ^^ ^^^ radio frequency in GHz, Fo{t) is the 
flux at 1 GHz in mjy and a{t) is the power-law index. 
The parameters Fo{t) and a{t) were characterized as a 
sum of Chebyshev polynomials of up to the 6-th power: 
Foit) = ^o,j Chebj, a{t) = ^bj Chebj. We note that 
such a characterization allows deriving the power-law index 
even if the data at different frequencies are not taken exactly 
at the same time. 

The fits to the radio LCs and the derived flux nor- 
malization (Fo) and power-law index (a) are shown in 
Fig. [71 Note that the variations of the power-law index are 
simila r to those derived in t h e original analysis of these 
data ("Mancheste r et all |2002| : Istavelev-Smith et all l2007l : 
[Zanardo et al. 2010|). Ne xt we inter pret these results by con- 
sidering the basic model (|Bellll 19781 ) for diffusive shock accel- 
eration of relativistic electrons and th e theory of synchrotron 
emission (jRvbicki fc Lightmanll 19791 ). 

In the framework of the DSA model, the resultant rel- 
ativistic particles have a power-law distribution in en ergy 
and t he corresponding power-law index is (eq. [12] in iBelJ 
[l97i): 



(2 + dj^) + dj^{2V^/Vsh - I/Ma) 
{dju - 1) - dju(V„/Kh + 1/Ma) 



(4) 



where dju is the density jump at the shock front, Vsh is 
the shock velocity, Vm is the velocity of the turbulent waves 
downstream and Ma is the Alfven Mach number of the 
shock. We note that for strong adiabatic shocks (as ex- 
pected in SNR 1987A) the density jump is dju = 4 for 
gas (plasma) with adiabatic index 7 = 5/3. Thus, three lim- 
iting cases are possible: (i) shocks with large Alfven Mach 
number and turbulent waves which have quickly isotropised 
(Vui <C Vsh, Ma ^ 1); (ii) strong shocks with similar turbu- 
lent waves as above but with Alfven Mach number not much 
larger than unity (dju = 4, V'^ < Vsh, Ma > 1); (hi) strong 
shocks with large Alfven Mach number and turbulent waves 
that have not isotropised yet {dju ~ 4, Vw/Vsh 7^ 0; Ma 3> 

The corresponding results for these three limiting cases 
are shown in the two upper panels of Figure [8] We made 
use of the relation between the spectral index of the 
resultant synchrotron emission from relativistic electrons 
with a power-law distribut ion in energy, a = —{p -— l)/2 
(|Rvbicki fc Lightmanlll979l ). In case (i), we see that the den- 
sity jump at the shock(s) in SNR 1987A is smaller than its 
canonical value (dju ~ 4) for strong shocks. This may in- 



dicate that efficient particle acceleration influenced the hy- 
drodynamic structure of the shock(s). If so, the derived val- 
ues for the density jump correspond to the subshock while 
the total compression (for the downstream versus upstream 
gas) is greater than 4 (e. g., for the case of SNR 1987A see 



iBerezhko fc Ksenofonto'^ (2OO6I : iDuffv et al.lll995h . In case 
(ii), we see that the Alfven Mach number has a tendency to 
increase with time. We note that the opposite is expected 
if the blast wave penetrates further on in the HII region 
near the inner ring and the magnetic fleld is subject to the 
frozen-in conditions (being proportional to the density). In 
case (iii), we see that after an initial increase the turbulent 
waves velocity is decreasing which may indicate that they 
evolve towards equilibrium — that is they tend to isotropize 
with time. 

On the other hand, the flts to the radio LCs of 
SNR 1987A tell us the time evolution of the normalization 
for its flux (_Fb) which in turn can be used to deduce some 
valuable information about the relati vistic electrons, mag- 
netic fi eld etc. From equation (6.36) in lRvbicki fc LightmarJ 
(Il979[ ) and for a distance of 50 kpc to SNR 1987A we can 
write: 



Fo = 4.272 X 10 



^^nZ>VB'-'^r^r2CSi-2a) 



2(1 -a) 



(5) 



,{'■) 



where _Fb is in mJy; n\J^ is the total number density of the 
relativistic electrons; V = V{t)/Vii = (r/_Riji„g)^A is the 
volume occupied by the shocked gas (and relativistic elec- 
trons), that is between the blast wave and the contact dis- 
continuity, in relative units and A = Ar/r is the relative 
thickness of that volume; B is the magnetic field strength 
in Gauss (B = -Bq-^ - frozen-in conditions, the subscript 
denotes parameters at the inner radius of the HII region); 
Co = (27r7nec/3e) x 10^ — 119.119 is a constant {me,e are 
the electron mass and charge, c is the speed of light, while 
the factor 10^ comes from using the radio frequencies in 
GHz); Fi = F(f + if) and F2 = r(f - ^) as F is the 
gamma function. 

Thus, from eq. ([SJ and the fit results for the radio LCs, 
we can derive the time evolution of the mean number den- 
sity of relativistic electrons. We recall that the radius of the 
blast wave, and the preshock density profile are those used 
in the fits to the X-ray CCD spectra of SNR 1987A. The 
relativistic number density can then be used to estimate the 
pressure of relativistic electrons in the radio-emitting region. 
Figure |8] presents results for a nominal set of values for the 
preshock nucleon density and magnetic field at the inner 
radius of the HII region (no = 100 cm~^ , Bo = 1 mG) 
and thickness (e.g., A « 0.1). We recall that the gas pres- 
sure behind the shock is pg = 3/4poVsh {po ~ l.Qmpno for 
the SNR 1987A abundances, rup is the proton mass). We 
see that neither the pressure of relativistic electrons nor the 
magnetic field pressure is a significant fraction of the ther- 
mal gas pressure. This indicates that the model of standard 
strong shocks can explain both the observed X-ray and ra- 
dio emission from SNR 1987A. Also, the small magnetic 
pressure is a sign of large Alfven Mach numbers, thus, case 
(iii) mentioned above seems a realistic option for the corre- 
sponding physical picture, and future observations will show 
us if this picture may evolve to being intermediate between 
cases (ii) and (iii). 

From eq. ([5]) and the fit results for the radio LCs, we 
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Figure 7. The radio LCs (pluses) of SNR 1987A (from lZanardo et al.ll201Cl ) and their model fits with eq.® (solid line). The model 
parameters, flux normalization at 1 GHz (-Fo) and spectral index (a), are shown in the bottom two panels. The dashed lines represent 
the corresponding Icr errors from the fit. The rms from the entire fit is 3.9 mjy while its value for the individual frequencies is 4.8 mjy 
(1.4 GHz), 5.2 mJy (2.4 GHz), 2.8 mjy (4.8 GHz) and 2.4 mjy (8.6 GHz). 



can also estimate the time evolution of the total number 
of relativistic electrons in the shock region and some re- 
lated quantities. Namely, the total number of particles is 
■nlfJ^V = N{t) = N{to) + N{t){t - to), where N{t) is the 
time-averaged production rate of relativistic electrons. Since 
the relativistic particles are accelerated at the shock front 
and then injected into the postshock region, we could also 



(r) 



(r) 



write: N{t) = N{to) + n\Ju2SAt, where niJ is the num- 
ber density of the injected relativistic electrons, S is the 
shock surface (5* oc [r / Rmng]'^) and accor ding to the DSA 
model: U2 = Kh/4 -I- Vn, for strong shocks (|Belllll978l ). Fig- 
ure [8] presents the time evolution of N{t) and n^„ based 
on the fits to the radio LCs and for case (iii). We see that 
after an initial steep increase of the newly born relativistic 
electrons, their average production rate follows a gradual 
decline. Note also that the number density of the newly in- 
jected relativistic electrons, normalized to the local preshock 
density, follows a similar pattern. We tentatively attribute 
this decline to decreasing efficiency of particle acceleration in 
the decelerating shocks. We believe that this behavior may 
be an important clue to the details of the particle accelera- 
tion mechanism that operates in SNR 1987A. Such details 
might be helpful to build a physically more sophisticated 
model of this system. 



6.2 Overall Consistency of the Global RSS 
Picture 

As a further constraint on this picture and as a check on its 
internal consistency, we can use the results from the fits 
to the CCD X-ray spectra of SNR 1987A. Namely, the 
preshock nucleon number density can be deduced from the 
emission measure of the blast wave and a density jump for 
strong adiabatic shocks {dju ~ 4) . For distance of 50 kpc to 
SNR 1987A the X-ray fits yield a nucleon number density at 
the inner radius of the HII region: no,x = (49.2 ± l)/\'sm9, 
cm~^, where d is the opening half-angle of the inner ring 
(6 = 7r/2 gives a complete spherical shell). Recall that we 
have a clear evidence based on the spatial-spectral effects in 
grating spectra that the X-ray emission in SNR 1987A is 
confined in the equatorial plane although the opening an- 
gle of the X-ray ern ission region is not yet well constrained 
(iZhekovet al.ll2005l . 2006; iDewev et aLll2008l : IZhekov et ai] 
2009J). On the other hand, the trip le-ring system observed 
in the optical (JBurrows et alJll995l ) is solid evidence for a 
highly asymmetric distribution of the circumstellar matter 
around the exploded star in the center of SNR 1987A. Al- 
though the origin of this asymmetry is not yet well under- 
stood, the scenario of colliding stellar winds that consid- 
ers a spherically symmetric BSG wind impinging onto an 
asy mmetric RSG is prob a bly th e most realistic model so 
far. iBlondin fc Lundgvistl ([1993) investigated this scenario 
in detail and their 'best' model, which correctly predicts 
various parameters of the circumstellar environment prior 
to the explosion of SN 1987A, requires that 50% of the RSG 
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Figure 8. The radio LC fit parameters of SNR 1987A. The Upper Left panel: the Alfven Mach number (sohd Une) and the density jump 
(dashed Une). The Upper Right paneh the turbulent velocity. The Lower Left panel: the time-averaged production rate of relativistic 
electrons (dashed line), the number density of the injected relativistic electrons (solid line) and the blast wave velocity (dotted line) as 
all these are normalized to their maximum value for clarity. The Lower Right panel: the ratios of the pressure of relativistic electrons 
(dashed line) and magnetic pressure (dotted line) to the thermal pressure behind the blast wave. See text for details. 



wind mass is coUimated within ~ 10° from the equatorial 
plane. We thus see that the nucleon density of the HII region 
that is deduced from the analysis of the CCD X-ray spectra 
is no,x = 118 ± 2.5, cm~^. This value is consistent with the 
one assumed for the estimates presented in Fig. [8] 

From these considerations we conclude that the case 
of strong adiabatic shock(s) is suitable for explaining both 
the X-ray and radio emission from SNR 1987A and that 
further analysis of the newly obtained data along the same 
lines might be indeed helpful in building a realistic picture 
of this object. It is worth noting, though, that this is in- 
deed the case provided the DSA of heavy particles is not 
as efficient as the one for the electrons in SNR 1987A. For 
example, if the relativistic protons gain the same amount 
of energy from the blast wave as supplied into relativistic 
electrons the picture will remain the same. But, if the DSA 
is very efficient and produces relativistic protons (or nuclear 
cosmic rays) and electrons with similar densities, the total 
pressure of all the relativistic particles will be comparable 
to the thermal gas pressure and the physical picture will be 
more similar to case (i) as described above. We recall that 
in such a case the blast wave in SNR 1987A will be affected 
by the feedback from the relativistic particles which will re- 
sult in changed physical parameters in the postshock region. 
A self-consistent treatment of these effects is needed which 
is beyond the scope of this study. We only note that the 
kinematics of the postshock (X-ray emitting) plasma will 
be influenced most, thus, a continuing monitoring of the X- 



ray emission from SNR 1987A with high spectral resolution 
(grating spectra) will be crucial in this respect. 

Some deceleration in the bulk gas velocities was found 
between 2004 and 2007 from the grating spectra with very 
good photon statistics. The most important kinematic result 
from these data, however, was the flnding that the bulk gas 
velocities inferred from the X-ray line profiles are too low 
to account for the p ostshock plasma tem pera tures derived 
from the spectral fits llZhekov et al.ll2005l . 2006: lDewev et alj 
'200^; IZhekov et all '2009)). On the other hand, if eflicient 
DSA operates in SNR 1987A, the total gas compression 
downstream from the shock front (i.e., in the X-ray emitting 
region) will be greater than its canonical value {dju,tot > 4). 
Such a result would imply that the bulk gas velocities in 
the space wifi be greater: Vbuik = y-'°'~ Vsh (e.g., 



also IChevalieilll983h . If so, the X-ray spectral lines would 
be broader, in confiict to observations. These considerations 
support our assumption that the DSA in SNR 1987A is not 
very efficient. 

This assumption gets additional support from the fact 
that the X-ray spectrum of SNR 1987A is dominated by 
emission from thermal plasma. Moreover, if, using the re- 
sults from the fits to radio LCs {Fo,a), we extrapolate the 
NT radio spectrum to X-rays (equivalent to assuming effi- 
cient particle acceleration), the observed flux in the (3 - 10 
keV) energy band is smaller than the one expected from rel- 
ativistic electrons for the SNR 1987A age larger than 5,500 
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DAE. This result confirms that the DSA is not very efficient 
in producing relativistic particles with very high energies. 

In the analysis presented here, we have argued that the 
physical picture in which the X-ray and NT radio emission 
from SNR 1987A originate in one and the same shock struc- 
ture is realistic and capable of describing the observations. 
But, there is one important detail that emerges from the 
analyses of the X-ray and radio images of this object and 
needs more attention. If the picture adopted here is correct, 
the LCs in X-rays and radio may differ but the angular size 
of SNR 1987A in these spectral domains must be about the 
same. However, the image analyses so far indicate that the 
X-ray size of SNR 1987A is consistently smaller than its ra- 
dio size. Because the object size is not measured directly but 
only after applying image- analysis techniques, such a differ- 
ence might be attributed to the different methods used in 
analyzing the X-ray and radio images (Gaensler ct al. 2007i). 
On the other hand, a recent and m ore elabor a te a nalysis 
of the radio images of SNR 1987A iNg et all (120081 ) find 
that for times beyond 6,000 days after the explosion (6,000 - 
7,600 DAE) the radius of the radio remnant is greater than 
0.83 arcsec, i.e., greater than the radius of the inner ring 
(see their Table 3). For the same period of t ime, the X-ray 
remn ant still lies within the equatorial ring (jRacusin et al.l 
I2OO9I ). 

Such an evolution of the radio size of SNR 1987A could 
be indeed important since it may indicate that the blast 
wave has propagated beyond the equatorial ring or the radio 
emission has some contribution from regions above and be- 
low it. The former was concluded by Ne et al. (2009) from 
their analysis of a single Chandra HRC (High Resolution 
Camera) image taken at 7736 DAE. They also derived a 
large r X-ray radius of SN R 1987A compared to the results 
from lRacusin et al.l (|2009l ). But if this were the case, the ra- 
dio (and or X-ray) emission should have started to decline 
similarly to the radio (X-ray) evolution soon after the ex- 
plosion when the supernova ejecta was interacting with the 
BSG stellar wind. This is to be expected in the framework 
of the colliding stellar wind scenario for the origin of the 
ring system in this object because the gas beyond the inner 
ring is that produced by the wind of the red supergiant pro- 
genitor of SNR 1987A. Therefore, the emission measure in 
radio and X-rays must decrease with time (preshock num- 
ber density will follow 1/r^ dependence). No indications of 
fiu x decline are seen both in the radio and X-ray LCs (e.; 
see iNg et al.ll2008l : IZanardo et aLll2010l : iRacusin et ai]l2009 
And in order to explore the possibility for contribution from 
shocks above and below the equatorial plane, we may need 
to change our currently used geometrical models. Namely, 
models that assume deviation from a simple geometry of a 
torus, equatorial belt etc. might be a more realistic presen- 
tation of the actual interaction of the blast wave with the 
inner ring. We believe that the improving quality of the X- 
ray and radio data allow for adopting such models which in 
turn may provide us with valuable details about the physics 
of the interaction phenomenon. 

An interesting aspect of our study is about what we can 
expect in the future, provided the physical picture adopted 
here is correct. As mentioned earlier (see ij 13.11 and Fig. [l|, 
the expansion velocity curve is successfully fitted for density 
profiles that may have difi'erent density contrast between the 
HII region and the maximum density of the smooth compo- 



nent at the radius of the inner ring. This results in a consid- 
erably different time {tmax) when the blast wave will reach 
the radius of the equatorial ring and the X-ray (as well as 
the radio) emission will be at its maximum. We fitted all 14 
CCD spectra simultaneously with the global RSS model for 
the cases with tmax between 10,000 and 15,000 DAE (with 
an increment of a 1000 DAE) . These fits had similar quality 
(X^ = 1696 - 1734) as our basic case {tma^ = 13, 000 DAE) 
discussed here in detail. As anticipated, it is not possible 
to discriminate between these cases just based on the avail- 
able data and only future studies of the expansion velocity 
curve and other X-ray properties will help us in this respect. 
But it is worth mentioning that all these fits converge in the 
following. 

Since the absolute value of the emission measure of 
the blast wave is constrained by the observed X-ray fiux, 
in all the cases the nucleon number density near the in- 
ner radius of the HII region had values that are about the 
same: no,x ~ 100 cm~^ (as before we assumed an open- 
ing half-angle of ~ 10°). From this it follows (see Fig [TJ 
that the maximum density of the smooth component will 



be URing 



W 



10 cm for the cases with t„ 



be- 



tween 10,000 and 15,000 DAE. On the other hand, in all 
these cases we need to have dense clumps with about the 
same typical number density {riciump ~ 10'* cm^'') in or- 
der to match the observed X-ray spectra. We note that the 
clump density (in fact its contrast to the smooth compo- 
nent) controls the properties of the transmitted shocks (low 
temperature plasma) . Interestingly, an indication of similar 
preshock number density (~ lO** cm~'^) was found also for 
the first hot spot discovered in the optic al and UV . Usini 
the line ratios of various forbidden lines, |Pun et al.l (|200! 
showed that the optical emission of Spot 1 comes from re- 
gions with densities of n^ ~ 10^ cm~^ downstream radia- 
tive shocks that compressed the preshock gas by a factor of 
^ 100. Thus, it seems conclusive that the HII region consists 
of dense clumps with similar properties and the density of 
its smooth component is not well constrained but it will be 
so depending on the time when the X-ray emission reaches 
its maximum. 

It is worth noting that although our global RSS model 
is successful in matching both the expansion velocity curve 
of SNR 1987A and the evolution of its X-ray emission, it 
neglects the role of the reverse shock for the latter. In fact, 
the emission from the entire two-slab region (see §[2]), that 
is of the shocked material between the forward and reverse 
shocks, is represented by the X-ray emission behind one 
fast shock: the evolving blast wave interacting with dense 
clumps. The validity of this seems justified given the good 
fits of our model to the X-ray spectra. So it could well be 
that the reverse shock produces unusually cool plasma, as it 
is a shock going through an almost neutral medium, in which 
magnetic turbulence, needed for coUisionless shock heating, 
is damped. However, there are enough uncertainties in our 
model, but also concerning the general knowledge about col- 
lisionless shocks to remain cautious. We only note that the 
physics of the reverse shock in SNR 1987A is quite com- 
plicated to be addressed properly in a simplified physical 
picture we have adopted in this work. 

As indicated by th e infrared observations (e.g. Fig. 18 
in lBouchet et al.ll2006l ) , the SN ejecta has become neutral 
already at times before the SNR 1987A reappearance in 
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the radio and X-rays (~ 1,000 DAE). On the other hand, 
dete ction of Ha and Lyg emission from th e reverse shock 
(e.g. iMichael et all l2003l : iHeng et al.l I2OO6I I is a clear sign 
for presence of both neutral and ionized species. Thus, the 
reverse shock in SNR 1987A is not a standard collisional 
or coUisionless shock in plasma but instead its shock front 
likely consists of a relatively broad transition zone where the 
ionization state changes from neutral (preshock ejecta ma- 
terial) into ionized (compressed ejecta plasma). As a result, 
it might well be that the electrons are not preheated in the 
shock front and their temperature is determined by the slow 
Coulomb collisions with protons and ions downstream the 
reverse shock which may not result in strong X-ray emission. 

But, there is an interesting transformation of the re- 
verse shoc k that we may witn ess in the near future. As sug- 
gested bv ISmith et al.l (120051 ). the rapidly increasing X-ray 
flux from SNR 1987A should reach a level beyond which it 
will be capable of ionizing the neutral gas before it crosses 
the reverse shock. When this happens, the broad Ha and 
Lya emission will vanish and most importantly the reverse 
shock will transform from its current state into a coUisionless 
shock, provided a magnetic field is present in the supernova 
ejecta. 

The reverse shock transformation will have an inter- 
estin g impact on the r a dio e mission from SNR 1987A as 
well. [Manchester et al.l (|2005r ) have suggested that the ra- 
dio emission originates from the region behind the reverse 
shock. Park et al. (2005, 2007) noticed the close similarity 
in the shape of the light curves in the hard X-rays and radio 
and suggested that they may have a common origin. But, no 
signs of non-thermal X-ray emission form SNR 1987A have 
been found so far and we argued here that the global RSS 
picture is capable of explaining the radio and X-ray prop- 
erties of SNR 1987A by assuming that the radio emission 
originates from the region behind the blast wave and not 
behind the reverse shock. On the other hand, once the su- 
pernova ejecta gas has been ionized before entering the re- 
verse shock, the shock will become coUisionless which will 
result in an efficient DSA of relativistic particles. After that 
moment, we can expect a more rapid increase of the non- 
thermal emission from SNR 1987A . 

We note that although this picture is only qualitative 
and thus bears many uncertainties if it is correct, two ob- 
servational evidences should be found in the future that can 
serve as a test for its validity. And it is important to note 
that these events must happen almost simultaneously: (i) the 
broad Ha and Lya emission from SNR 1987A will vanish; 
(ii) there will be a turn-up in the non-thermal radio emission 
from this object. 

6.3 SNR 1987A vs. Young Supernova Remnants 

Many young supernova remnants (Gas A, Kepler, Tycho, 
SN 1006, RCW 86) have been observed to contain non- 
thermal X-ray emission (e.g. lAllen et al.lll999l ). These young 
remnants share other common features: non-thermal radio 
emission, signs of efficient (nuclear) cosmic ray (CR) accel- 
eration and as a result from that they likely have ampli- 
fied magnetic fields (for a discussion see Volk et al. 2005; 
Vink 2008) We note that the key signature of efficient DSA 
of relativistic particles in these objects is the presence of 
non-therm,al X-ray emission. As already mentioned (^ 16. 2|) . 



SNR 1987A does not have that (or it is very weak) and its 
X-ray spectrum is successfully represented by thermal emis- 
sion from shock heat ed plasma as analyses of both Chandra 
(:Michael et al.ll2002l : Park et al. 2002, 2004, 2006; Zhekoy 
et al. 2006 200 9) and XMM-Newton (JHaberl et al.l I2OO6I : 
iHeng et al.ll2008r ) spectra have shown. This is also confirmed 
by our current analysis of the evolution of the Chandra CCD 
spectra of SNR 1987A . 

It is then suggestive that the lack of non-thermal X-ray 
emission from SNR 1987A could be a signature of inefficient 
DSA of relativistic particles (protons and nuclear CR) and 
it might simply be explained by the fact that SNR 1987A is 
very young. Actually, its age is more than one order of mag- 
nitude smaller than that of the young SNRs showing those 
common features mentioned above. We can speculate that 
some time is probably needed before a SNR reaches a level 
when the DSA becomes very efficient and SNR 1987A has 
not yet entered that stage. On the other hand, old SNRs 
are not expected to have strong non-thermal X-ra y emission 
since their shock velocities are deccelerating (e.g. iVolk et al.l 
[2OO5I). This can be an alternative explanation for the case 
of SNR 1987A that shows c lear signs of shock d eccelera- 
tion at times > 5,500 DAE (jRacusin et al.ll2009l : see also 
Fig. [T] and the lower right panel in Fig. [3| . So, due to the 
presence of the inner equatorial ring and its HH region (a 
CSM with increasing outward density), the very young ob- 
ject, SNR 1987A , might mimic the evolution of old SNRs. 
But in either case, the fact of missing (or very weak) non- 
thermal X-ray emission probably indicates inefficient DSA 
in SNR 1987A . 

Thus, the assumption of not very efficient DSA was 
adopted in our analysis here which also means that the mag- 
netic field strength cannot be derived in a self-consistent 
manner as can be done for young SNRs with non-thermal 
X-ray emission (e.g. Volk et al. 2005). Due to this, our es- 
timates (see § 16.11 also the lower panels in Fig. [SJ use some 
fiducial value for the magnetic field strength at the inner ra- 
dius of the HH region and we adopt the frozen-in conditions. 
The latter means that the field strength increases deeper in 
the inner ring. We note that in the interacting wind scenario, 
the origin of the magnetic field is likely related to the global 
field of the progenitor star during the BSG (and RSG) phase 
of its evolution. The adopted here value (Bo — 1 mC) sug- 
gests that not a strong surface field is required if the global 
structure of the stellar magnetic field was according to the 
Parker's model (|Parkeij[l958 '). 

We note that such not strong magnetic fields are likely 
present in massive hot stars and the clear sign for this is the 
non-th ermal radio emission detected from wide WR-(-0 bi- 
naries IjDoughertv fc Williamsll2000l ). This emission is iden- 
tified with the colliding wind shock region in these objects 
and WR 140 and WR 147 are two classical examples of those. 
Interestingly, there are no signs of non-thermal component 
in the X-ray spectra of these objects which are successfully 
matched by thermal emission from plas ma behind strong 
adiabatic shocks (e . g. for WR 140 see IZh ekov & SkinneiJ 
[2OOOI: iPoUock et all I2OO5I : for WR 147 see Zhekoy, :20o3; 
IZhekov fc Parkll2010r ). Thus, a physical picture where strong 
adiabatic shocks produce thermal X-rays and the corre- 
sponding DSA of relativistic particles is not very efficient 
does not seem to be anything unique and it might well be 
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a reasonable assumption for the current evolution stage of 
SNR 1987A . 

Finally, the ongoing monitoring of the X-ray emission 
from SNR 1987A will be very helpful in further constrain- 
ing the physical picture described here. First, the continu- 
ous brightening of SNR 1987A will allows us to establish 
with higher certainty the level of non-thermal X-ray emis- 
sion from this very young supernova remnant. We note that 
if strong (and increasing) non-thermal emission is detected, 
it will require considerable changes in the physical picture as 
we described it in this work. Second, the data from the mon- 
itoring observations will also provide us with detailed spec- 
tral information (grating data) that will reveal the evolution 
of the kinematics of the X-ray emitting plasma. We note that 
the RSS picture was proposed to explain the relatively nar- 
row spectral lines in the grating spectra of SNR 1987A. As 
seen from Figs. [5] and (6] the contribution from the reflected 
shock to the total observed fluxes is not dominant even at 
times after 6,000 DAE which might be a deficiency from 
the poor spectral resolution in CCD spectra. Upgrading the 
available CCD data base with grating spectra, having good 
quality and spanning over long enough period of time, will 
allow us to test the global RSS model more rigorously by 
confronting its predictions with the actual evolution of the 
gas kinematics of the X-ray emitting plasma. This will also 
justify amendments to the model by taking into account the 
X-ray emission from the region behind the reverse shock. In 
turn, we can further improve our understanding of the excit- 
ing phenomenon of the birth of supernova remnant 1987A. 



ray LC reported bv IPark et al.l (|2005f ). Also, the emission 
measure of the hot plasma in the transmitted shock (as well 
as in the reflected shock) is increasing steeper than that 
for the blast wave for times after ~ 6400 DAE. All this is 
indicative of an ongoing penetration of the blast wave deeper 
in the equatorial ring debris. 

4. Results from the X-ray analysis were used to explore 
the possibility that the X-rays and the non-thermal radio 
emission from SNR 1987A originate from the same shock 
structure (the blast wave). Using a simplified analysis, it 
was shown that this can indeed be the case, provided the 
particle acceleration mechanisms, that operate in the shocks 
in SNR 1987A, are not very efficient in producing relativistic 
particles with very high energies. 

5. As a next step, the global RSS model could be used 
to obtain further constraints on the reflected-shock struc- 
ture picture in SNR 1987A by analyzing X-ray spectra with 
good spectral resolution once the available data cover a long 
enough time interval of the SNR 1987A evolution. We be- 
lieve that this will allow us to improve our understanding of 
the underlying physics and to build a more realistic picture 
and model of this fascinating phenomenon: the birth and 
evolution of supernova remnant 1987A. 
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7 CONCLUSIONS 

In this study, we continued to explore the validity of the 
reflected-shock structure picture in SNR 1987A. This pic- 
ture assumes that the blast wave is interacting with the 
HII region in vicinity of the inner ring and when it encoun- 
ters dense clump(s) a shock is transmitted into the clump(s) 
and a reflected shock forms that additionally compresses the 
shocked gas behind the blast wave. We adopted an improved 
version of the reflected-shock model (global RSS model) to 
test the RSS picture in SNR 1987A by making use of CCD 
spectra from the monitoring program with Chandra. The 
main results from this analysis are as follows. 

1. The global RSS model is ca pable of matching b oth 
the X-ray expansion velocity curve l|Racusin et al.ll2009l ') as 
well as the X-ray spectra of SNR 1987A over its evolution 
for the last ten years or so. Moreover, extrapolating back in 
time the X-ray flux, the model is able to match the observed 
ROSAT fluxes from the time of the supernova reappearance 
in X-rays (1000-3000 DAE). 

2. The evolution of the mean electron temperature be- 
hind the shocks shows a decrease with time of the electron 
temperature in the blast wave, while a slight increase with 
time is noticed for this parameter in the transmitted shock. 
The mean ratio of the electron to the mean plasma tempera- 
ture is required to increase with the decrease of the velocity 
of the blast wave. 

3. At early times, the soft (0.5-2 keV) X-ray flux of 
SNR 1987A was dominated by the blast wave but for times 
after ~ 6000 DAE the transmitted shock became dominant. 
This transition coincides with the steep upturn in the X- 
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